Kinetic theory of electromagnetic ion waves in relativistic plasmas 
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y-H , A kinetic theory for electromagnetic ion waves in a cold relativistic plasma is 

derived. The kinetic equation for the broadband electromagnetic ion waves is coupled 
co ■ to the slow density response via an acoustic equation driven by ponderomotive force 

ci 

like term linear in the electromagnetic field amplitude. The modulational instability 
growth rate is derived for an arbitrary spectrum of waves. The monochromatic and 
random phase cases are studied. 
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The interaction of intense radiation with plasmas is of fundamental importance in a wide 



variety of applications, such as inertial confinement fusion [1|, |2| and pulsar emissions [3] , and 

\o ■ 

£f) ' such interactions can lead to a range of instabilities, e.g. Brillouin and Raman scattering as 



00 ■ well as modulational instabilities m, y, IfjLJaL In some of the instabilities, relativistic effects 
play an important role . In Refs. [lfj, |Tj| it has been shown that a novel type of ion 

O ' 

electromagnetic wave may propagate through a cold magnetized relativistic plasma. Such 

o ; 

waves may be of relevance in high-energy astrophysical environments, e.g. pulsars, as well 
as in intense laser interactions with high density targets. However, the coherence length 



of intense electromagnetic (EM) waves could often be short. Moreover, it is well-known 



that effects of partial coherence may be used to stabilize the propagation of EM pulses in 
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nonlinear dispersive media 



12l Q|. Thus, it is of interest to analyze the effects of partial 
coherence of EM ion wave propagation. 

In this paper, the effects of partial coherence ion EM waves in a magnetized relativistic 
plasma is analyzed. In particular, a kinetic equation describing the evolution of the EM 
quasi-particles is derived and the general nonlinear dispersion relation is obtained. We 
study the effects of partial coherence using a random phase background of EM ion waves, 
leading to a Lorentz distribution for the quasi-particles. It is shown that the spatial spectral 
broadening induced by the random phase gives rise to a reduction of the modulational 
instability growth rate, in addition to the appearance of a new short wavelength instability 
region. 

Let us consider large amplitude circularly polarized electromagnetic ion waves, with the 
electric field amplitude E , propagating parallel to a constant magnetic field B = B z in a 
cold electron-ion are situations in 

which the ions are strongly relativistic. For such electromagnetic waves, the linear dispersion 
relation is of the form 

k 2 r 2 UJ 2 - 

1«^±-^, (1) 

where k is the wavenumber, uj is the wave frequency, cu pi = (e 2 n 0i /eQm i ) 1 ^ 2 is the ion plasma 
frequency, uje% = eE /cmi, e is the magnitude of the electron charge, n 0i is the unperturbed 
ion density, eo is the vacuum dielectric constant, rrii is the ion rest mass, and c is the speed of 
light in vacuum. The + (— ) in ((TJ) correspond to \pi\/rriiC > u C i/u (\pi\/m,iC < u} ci /u), where 
Pi is the ion particle momentum and u C i = eB /mi is the ion cyclotron frequency. If the 
plus sign is used in the dispersion relation we obtain the so called second ion-cyclotron 
wave, whereas the minus sign gives the ion helicon wave. 

The nonlinear evolution of electromagnetic ion waves, with amplitude E, in a cold quasi- 
neutral relativistic plasma is given by 

d d\ „ v'd 2 E f\E\-E 5n\ „ 
where the + (— ) refers to the + (— ) in the linear dispersion relation (0), and the slow density 
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response is obtained from 



d 2 2 d 2 \6n d 2 [ \E\-Eq . 

v sTOT — = ± fe F • ( 3 ) 



at 2 W J n ^ 2 V 
Here v g = dui/dk is the group velocity, v' = d 2 u/dk 2 is the group velocity dispersion, 
Sn is the ion density perturbation, no is the constant background density, v s = (T e /rrii) l l 2 
is the ion sound speed, and (3 = (c 2 /2)(E /cB ). Thus, the above system of equations 
is effectively a quadratic nonlinear Schrodinger equation, in contrast to the regular cubic 
nonlinear Schrodinger equation. 

In the quasi-stationary limit, we may integrate Eq. (jSJ) to obtain the quadratic nonlinear 
Schrodinger equation 



; dE v' g d 2 E Y i , j3\ \E\ 



! 97 + W ± H 1± U^ = ' (4) 

where we have transformed to a co-moving frame t = t and ( = z — v g t. 

With the ansatz E = (Eq + 5E) exp(i<p), where 5E -C Eq and <C 1, we linearize Eqs. (J2J) 
and (JHJ). For a harmonic dependence, i.e. 5E, <fi oc exp(iKz — ittt), we obtain the dispersion 
relation 

(n 2 - v 2 K 2 ) [(fl - t^) 2 ± ±u;^ 2 - i<i, 4 ] = \^v' g K\ (5) 

consistent with the results in Refs. jll| and jl^j]. In the quasi-stationary limit (cf. Eq. (J3J)), 
we can solve for Q to obtain the growth rate T = lm(Q) according to 

1/2 

(6) 

For positive group velocity dispersion, we see that the second order dispersive term competes 
with the nonlinear term, giving the characteristic modulational instability growth rate curve 
(cf. Figs. 1 and 2). 

We are now interested in analyzing the effects of partial coherence of the electromagnetic 
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ion waves. For this purpose, we introduce the Wigner function 
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p(t,z,p) = ±- J d(e^E*(t,z + (/2)E(t,z-(/2), (7) 
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such that 

\E\= Udpp(t,z,p)J . (8) 

We note that the Wigner method [19], as well as the equivalent mutual coherence 
method |2(3j], have been used to analyze the modulational instability of the cubic nonlinear 
Schrodinger equation relevant for e.g. nonlinear optics. Moreover, in plasma applications 



the Zakharov equations have been analyzed using the above method 2l( in order to obtain 
the statistical dynamics and Landau like damping of Langmuir waves. 

Applying the time derivative to the definition Q and using (J2]L we obtain the kinetic 
equation 



dp 
dt 



i , \ dp , n (\E\ 5n\ . 1 d d \ 



where the sin-operator is determined in terms of its Taylor expansion and the arrows denotes 
direction of operation. The kinetic equation Q, together with © and (jSJ), determines the 
evolution of broad band electromagnetic ion waves in cold relativistic plasmas. 

In order to analyze the modulational instability of the system (J3J), (JEJ), arid (J5jl. we let 
p(t, z,p) = po(p) + <5p exp(ii^2; — ifii) and 5n oc exp(ifT,2 — iQt), where 5p <C po- Linearizing 
Eqs. (J3J), (JHJ), and Q, we then obtain the nonlinear dispersion relation 



^-^{^ ?w^w ) J dp « - (t + ' (10) 

where £o = (J dppo) 1 ^ 2 . This dispersion relation generalizes (j3J) to the case of arbitrary 
spatial spectral background distributions po- 

In the case of a monochromatic spectral distribution, po(p) = Eq5(p), we retrieve the 
dispersion relation ( |5|) . as expected. However, if the background distribution po has a finite 
spectral width, the dispersion relation is altered. Next, we look at the case of a random 
phased background electromagnetic ion wave. This will give rise to the spectral distribution 
in the form of the Lorentzian 



5 




(a) (b) 

FIG. 1: The growth rate T = Im(fi) plotted as a function of the normalized wavenumber K for 
the parameter values v = 0.33, = 1, and v' = 0.5, with + in Eq. (J2J). In (a) we have A = 0, while 
(b) has A = 0.2. The reduction in growth rate can be seen from (a) to (b). 

where A is the spectral width of the distribution. With this, the nonlinear dispersion relation 
(fTUj) becomes 

(fi 2 - v 2 s K 2 ) [[fl - (v g - i v' g A)Kf ± \uv> g K* - \vfK*] = \u;(3v' g K\ (12) 

As A — > 0, we obtain from above the monochromatic dispersion relation (j5J). Moreover, the 
effect of the spectral broadening is to introduce a damping of the perturbation modes. 

Next, we analyze the modulational instability properties of the dispersion relation (|12j). 
We introduce the dimensionless variables Q = VL/uj, K = v g K/u, A = v g A/u, v = v s /v g , 
v' = uv' g /Vg, and /3 — /3/v g . In Figs. 1 (+ in Eq. Q) and 2 (— in Eq. ©) we have displayed 
the growth rate T = lm(Q) for some typical parameter values. Using v = 0.33, (3 = 1, and 
v' = 0.5, it can be seen that a finite spectral width A gives rise to a reduced growth rate. 
On the other hand, the ^-region, where the instability occurs, is enlarged. We also note 
that the +-mode has a larger growth rate as compared to the — mode in Eq. @. 

To summarize, we have analyzed the effects of a partial coherence of circularly polarized 
electromagnetic ion waves in relativistic plasmas. In particular, the modulational instability 
growth rate was found for both coherent waves and for waves with a random phase. It was 
shown that the effect of partial coherence is to stabilize the propagation of the aforemen- 
tioned waves, while broadening the possible instability wavenumber region. The latter could 
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(a) (b) 

FIG. 2: The growth rate T = Im(fi) plotted as a function of the normalized wavenumber K for 
the parameter values v = 0.33, = 1, and v' = 0.5, with — in Eq. l[2*|). In (a) we have A = 0, while 
(b) has A = 0.2. The reduction in growth rate can be seen from (a) to (b). 

lead to the formation of short wavelength nonlinear structures due to partial coherence. 



Acknowledgments 

This research was partially supported by the Swedish Research Council. We thank L. 
Stenflo for helpful discussions and valuable insights. 



[1] G. A. Mourou, T. Tajima, and S. V. Bulanov, Rev. Mod. Phys. 78, 309 (2006). 
[2] M. Marklund and P. K. Shukla, Rev. Mod. Phys. 78, 591 (2006). 

[3] V. S. Beskin, A. V. Gurevich, and Ya. N. Istomin, Physics of the Pulsar Magnetosphere 
(Cambridge, 1993). 

[4] M. Y. Yu, K. H. Spatschek, and P. K. Shukla, Zh. Naturforsch. A 29, 1736 (1974). 

[5] P. K. Shukla, M. Y. Yu and K.H. Spatschek, Phys. Fluids 18, 265 (1975). 

[6] R. P. Sharma and P. K. Shukla, Phys. Fluids 26, 87 (1983); P. K. Shukla and L. Stenflo, ibid. 

28, 1576 (1985); G. Murtaza and P. K. Shukla, J. Plasma Phys. 31, 423 (1984); P. K. Shukla 

and L. Stenflo, Phys. Rev. A 30, 2110 (1984). 
[7] P. K. Shukla, N. N. Rao, M. Y. Yu, and N. L. Tsintsadze, Phys. Rep. 135, 1 (1986). 



7 



[8] N. L. Tsintsadze and L. Stenflo, Phys. Lett. A 48, 399 (1974). 

[9] C. E. Max, J. Arons, and A. B. Langdon, Phys. Rev. Lett. 33, 209 (1974). 
[10] N. Y. Kotsarenko, Soviet J. Plasma Phys. 3, 197 (1977). 
[11] L. Stenflo and N. L. Tsintsadze, Astrophys. Space Sci. 64, 513 (1979). 

[12] Y. Kato, K. Mima, N. Miyanaga, S. Arinaga, Y. Kitagawa, M. Nakatsuka, and C. Yamanaka, 

Phys. Rev. Lett. 53, 1057 (1984). 
[13] M. Koenig, B. Faral, J. M. Boudenne, D. Batani, A. Benuzzi, and S. Bossi, Phys. Rev. E 50, 

R3314 (1994). 

[14] N. L. Tsintsadze, N. A. Papuashvili, E. C. Tsikarishvili, and L. Stenflo, Physica Scripta 21, 
183 (1980). 

[15] E. P. Wigner, Phys. Rev. 40, 749 (1932); 

[16] J. E. Moyal, Proc. Cambridge Philos. Soc. 45, 99 (1949). 

[17] Yu. L. Klimontovich, The Statistical Theory of Non-Equilibrium Processes in a Plasma (Perg- 

amon Press, Oxford, 1967). 
[18] J. T. Mendonga, Theory of Photon Acceleration (IOP Publishing, Bristol, 2001). 
[19] D. Anderson, B. Hall, M. Lisak and M. Marklund, Phys. Rev. E 65, 046417 (2002); B. Hall, 

M. Lisak, D. Anderson, R. Fedele, and V. E. Semenov, ibid. 65, 035602 (2002). 
[20] M. Soljacic, M. Segev, T. Coskun, D. N. Christodoulides, and A. Vishwanath, Phys. Rev. 

Lett. 84, 467 (2000). 

[21] R. Fedele, P.K. Shukla, M. Onorato, D. Anderson, and M. Lisak, Phys. Lett. A 303, 61 (2002). 



